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Section  1 


INTRODUCTION 

This  is  the  third  of  a series  of  reports  on  the  development 
of  finite-element  analysis  procedures  for  computation  of  linear 
elastic  stress  intensity  factors  associated  with  cracks  in  common 
aircraft  structural  details.  The  first  report  in  the  series  [1] 
reviewed  the  foundations  of  the  assumed-stress  hybrid  finite-element 
method  and  its  application  to  the  formulation  of  a special  element 
which  contains  a crack-tip  singularity,  but  which  is  also  compatible 
with  conventional  assumed-displacement  elements  along  its  boundaries. 
The  second  report  [2]  documented  some  additions  to  the  ASRL  FEABL-2 
software  for  substructure  analysis,  the  development  of  a new  hybrid 
element  to  connect  the  vicinity  of  a fastener  hole  to  a cartesian 
panel  geometry,  and  the  program  codes  for  several  pieces  of  sub- 
structure which  were  put  together  to  model  a tension  panel  with  a 
single  fastener  hole.  Performance  tests  which  were  conducted  on 
the  new  hybrid  element  demonstrated,  among  other  things,  that  the 
element  could  function  accurately  at  a fastener  hole  subjected  to 
a cosine  bearing  pressure  distribution  applied  either  directly  to 
the  hybrid  element  or  to  the  inner  boundary  of  some  rings  of 
quadrilateral  elements  inserted  between  the  hybrid  element  and  the 
fastener  hole  surface.  In  the  present  report,  the  pieces  and  parts 
are  shuffled  and  presto! --two  new  finite-element  models  appear. 

They  are  tension  panels  with  either  a single  row  of  fastener  holes, 
or  a double  row  of  holes  on  staggered  centerlines.  A new  procedure 
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is  introduced  for  each  model:  SBL  and  DBL,  which  generate  the 

far-field  substructures.  However,  the  new  procedures  are  execu- 
tives which  merely  use  the  previously  developed  basic  substructures 
to  create  different  design  details  and  different  loading.  The 
fastener  hole  rows  are  now  located  along  the  top  edge  of  the  panel, 
and  tension  is  introduced  by  a cosine  bearing  pressure  distribution 
at  each  hole. 
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Section  2 


SR0V7  AND  DROW  PROGRAMS 

Parametric  codes  designated  SROW  and  DROW  have  been  developed 
to  generate  and  analyze  the  single-  and  double-fastener-row  panels. 
Programs  SROW  and  DROW  are  executive  controllers  which  direct  pro- 
cedures SBL  and  DBL /respectively , as  well  as  other  previously 
developed  procedures. 

2 . 1 Model  and  Input  Conventions  for  Program  SROW 

Fig.  1 illustrates  the  structure  which  is  modeled  by  program 
SROW:  a skin  panel  with  a horizontal  row  of  equally  spaced  fastener 

holes  near  the  top  edge.  The  centerline  of  the  row  is  placed  below 
the  edge  by  half  the  center-to-center  spacing  of  the  fasteners.  The 
panel  may  have  up  to  10  holes.  The  bottom  edge  of  the  panel  is 
restrained.  Load  is  introduced  by  fastener  bearing  forces  applied 
at  the  holes.  These  forces  are  represented  by  cosine  pressure  dis- 
tributions. The  bearing  forces  may  be  different  for  different 
fastener  holes.  The  model  assumes  that  one  fastener  hole  is 
damaged.  The  user  specifies  the  location  of  the  damaged  hole, 
whether  it  has  a single  crack  or  two  cracks  180  degrees  apart,  and 
the  range  of  angular  crack  positions  to  be  analyzed. 

Fig.  2 summarizes  the  input  data  card  conventions.  Five  or 
six  input  data  cards  are  required  to  supply  the  following  information: 
1 . Panel  Parameters 

NHOLES  = Number  of  fastener  holes  in  panel 
CLSPC  = Center-to-center  hole  spacing 
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RI  = Hole  radius 

IDMGD  = Position  of  the  damaged  hole 

LENGTH  = Total  length  of  the  panel 

THK  = Panel  thickness 

2 . Load  Parameters 

BFORCE.  = Total  bearing  force  to  be  applied  to  the 
1 ith  hole  in  the  row 

(Note:  these  data  are  continued  on  a second  card 

for  panels  with  NHOLES  > 8.) 

3 . Crack  Parameters 

A ( 1 ) = Length  of  first  crack 

A (2)  = Length  of  second  crack 

IP0S(1)  = Initial  position  number  of  first  crack 

IPOS (2)  = Final  position  number  of  first  crack 

4 . Material  Properties 

E = Young ' s modulus 

v = Poisson's  ratio 

5 . Print  Control  Parameters 

KT1  = Control  for  optional  FEABL-2  output 

KT2  = Control  for  procedure  SBL  output 

KT3  = Control  for  procedure  RING  output 

Any  value  from  1 to  10  may  be  chosen  for  NHOLES.  The  fastener 

hole  radius  is  restricted  to:  CLSPC/12  <_  RI  CLSPC/8.  The  fastener 

holes  are  considered  to  be  in  the  sequence  1,2,...,  NHOLES  from 
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left  to  right.  The  damaged  hole  is  defined  by  its  sequence  number: 

1 < IDMGD  £ NHOLES  (1) 

The  panel  width  is  computed  internally  as  NHOLES  * CLSPC  by  program 
SR0V7. 

Bearing  forces  for  each  fastener  hole  are  specified  individually 
on  the  second  data  card,  with  a continuation  card  required  if  the 
panel  has  9 or  10  holes.  The  bearing  force  values  may  be  zero  or 
positive.  Negative  values  (i.e.,  compressive  load  transfer)  can 
be  processed  by  the  program,  but  will  be  incorrectly  applied  as  a 
distributed  tension  on  the  upper  half  of  the  hole  rather  than  as 
a bearing  pressure  on  the  lower  half.  Therefore,  the  use  of  nega- 
tive bearing  forces  is  not  recommended  for  the  present  version  of 
program  SROW.  The  bearing  forces  are  converted  internally  to  con- 
sistent nodal  forces  which  model  a cosine  pressure  distribution 
statically  equivalent  to  the  bearing  force. 

Input  parameters  A(l),  A(2)  are  the  actual  crack  lengths,  as 
measured  from  the  edge  of  the  fastener  hole  to  the  crack  tip.  The 
cracks  are  placed  so  as  to  emanate  radially  from  the  fastener  hole. 

If  a case  with  only  one  crack  is  to  be  analyzed,  specify  A (2)  = 0.0. 
Crack  lengths  in  the  range: 

0 £ A(l) , A (2 ) <1.27 (RI)  (2) 

are  permitted.  The  initial  and  final  angular  positions  of  the 
first  crack  are  controlled  by  the  crack  position  numbers,  in 
accordance  with: 

6 initial  * «K)S(1)-1)A0 

® final  = (IPOSU)-DSe  (3) 

A0  = it/12  rad.  = 15  deg. 
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If  a second  crack  exists,  it  is  always  positioned  180  degrees  from 
the  first  crack.  Permissible  values  of  the  position  numbers  are: 


1 < IPOS  (1)  < n 

— — max 

IPOS(l)<  IPOS (2)  < n 

— — max 

where 


(4) 


n - 12 

max 

if 

A (2) 

O 

O 

"Ik 

n =24 

max 

if 

A (2) 

= 0.0 

Solutions  are  executed  automatically  for  the  angular  positions 

from  0.  . , . . to  0^.  in  15-deqree  increments.  A single  solution 

initial  final 

is  executed  if  IPOS (2)  = IPOS(l). 

The  print  control  parameters  may  be  used  to  delete  nonessential 
output  by  setting: 

KTi  = FORTRAN  unit  number  for  line 
printer  at  user's  computing 
facility . 

The  FORTRAN  unit  number  is  the  value  specified  for  the  line  printer 
in  a FORTRAN- IV  print  instruction;  e.g.: 

WRITE  (6,1000)  A,  B,  C 
(Unit  number  = 6) 

Any  other  values  for  KTi  (i  = 1,2,3)  will  permit  the  nonessential 
output  to  be  printed.  Full  output  is  recommended  for  the  initial 
trial  run  when  program  SROW  is  first  imported  to  a computing  facility. 
Deletion  of  all  three  options  is  recommended  for  production  runs. 


2 . 2 Model  and  Input  Conventions  for  Program  PROW 

Fig.  3 illustrates  the  structure  which  is  modeled  by  program 
DROW:  a skin  panel  with  two  horizontal  rows  of  equally  spaced 
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fastener  holes  near  the  top  edge.  The  upper  and  lower  rows  are 
designated  row  1 and  row  2,  respectively.  The  centerline  of  row  1 
is  placed  below  the  top  edge  by  half  the  fastener  center-to-center 
spacing.  Row  2 is  staggered  to  the  right  of  row  1 by  half  the 
center-to-center  spacing.  The  panel  may  have  up  to  8 fastener 
holes  per  row.  The  bottom  edge  of  the  panel  is  restrained.  Load 
is  introduced  by  fastener  bearing  forces  applied  at  the  holes. 
These  forces  are  represented  by  cosine  bearing  pressure  distribu- 
tions. A different  bearing  force  may  be  specified  for  each  hole. 
The  model  assumes  that  one  fastener  hole  is  damaged.  The  user 
specifies  the  location  of  the  damaged  hole,  whether  it  has  a 
single  crack  or  two  cracks  180  degrees  apart,  and  the  range  of 
angular  crack  positions  to  be  analyzed. 

Fig.  4 summarizes  the  input  data  card  conventions.  Six  input 
data  cards  are  required  to  supply  the  following  information: 

1 . Panel  Parameters 

NHOLES  = Total  number  of  holes  per  row 
CLSPC  = Center-to-center  fastener  spacing 
RI  = Hole  radius 

IROW  = Row  which  contains  damaged  hole 

IDMGD  = Position  of  damaged  hole  within  row 
LENGTH  = Total  length  of  panel 

THK  = Panel  thickness 

2 . Load  Parameters  (Row  1) 

BFORCE.=  Total  bearing  force  to  be  applied 
1 to  the  ith  hole  in  the  row 
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3 . Load  Parameters  (Row  2) 

BFORCE . = Total  bearing  force  to  be  applied 
1 to  the  ith  hole  in  the  row 

4 . Crack  Parameters 

A(l)  = Length  of  first  crack 
A (2)  = Length  of  second  crack 

IP0S(1)=  Initial  position  number  of  first  crack 
IP0S(2)=  Final  position  number  of  first  crack 

5 . Material  Properties 

E = Young's  modulus 

v = Poisson's  ratio 

6 . Print  Control  Parameters 

KT1  = Control  for  optional  FEABL-2  output 
KT2  = Control  for  procedure  DBL  output 

KT3  = Control  for  procedure  RING  output 

Any  value  from  1 to  8 may  be  chosen  for  NHOLES.  The  fastener 
hole  radius  is  restricted  to: 

CLSPC/12  < RI  < CLSPC/8  (5) 

The  fastener  holes  are  considered  to  be  in  the  sequence  1,2,..., 
NHOLES  from  left  to  right  in  each  row.  The  location  of  the  damaged 
hole  is  specified  by  row  number  IROW  = 1 or  2 and  position  number 
within  the  row  (Eq.  1) . The  panel  width  is  computed  internally  by 
program  DROW  as  CLSPC  * (NHOLES  + 0.5). 

Bearing  forces  for  each  fastener  hole  are  specified  individually 
on  the  second  data  card  for  row  1 and  on  the  third  data  card  for 
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row  2.  The  bearing  force  values  may  be  zero  or  positive.  Nega- 
tive values  (i.e.,  compressive  load  transfer)  can  be  processed  by 
the  program,  but  will  be  incorrectly  applied  as  a distributed  ten- 
sion on  the  upper  half  of  the  hole  rather  than  as  a bearing  pres- 
sure on  the  lower  half.  Therefore,  the  use  of  negative  bearing 
forces  is  not  recommended  for  the  present  version  of  program  DROW. 
The  bearing  forces  are  converted  internally  to  consistent  nodal 
forces  which  model  a cosine  pressure  distribution  statically 
equivalent  to  the  bearing  force. 

The  crack  parameters  and  print  control  parameters  are  defined 
in  the  same  manner  and  subject  to  the  same  restrictions  as  for 
program  SROW.  See  the  latter  part  of  Subsection  2.1  for  details. 

2 . 3 Required  Subprograms  and  Other  Features 

Programs  SROW  and  DROW  require  the  following  subprograms  for 
execution : 

1.  ASRL  FEABL-2  subroutines  ASMLTV , ASMSUB,  BOON,  FACT, 

ORK,  QBACK,  SETUP,  SIMULQ,  STACON,  and  XTRACT . 

2.  IBM  Scientific  Subroutine  Package  subroutines  MFSD 
and  SINV. 

3.  ASRL  elements  PCRK59  and  QUAD 4 . 

4.  Procedures  HOLEL,  LUG,  and  RING,  together  with  their 
associated  subroutines. 

5.  Procedure  SBL  (for  program  SROW)  or  procedure  DBL 
(for  program  DROW) . 

No  external  (disk  or  tape)  files  are  required  by  either  program. 
The  core  memory  requirements  are  as  follows: 

Program  SROW  Program  DROW 

344 

86 , 000. n 
247 , 760g 


272 

68 ,000. n 
204 , 640g 
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KBYTES 

Decimal  Words 
Octal  Words 


Execution  times  range  from  0.25  to  1.5  CPU  minutes,  depending  upon 
the  range  of  crack  parameters  investigated.  The  above  statistics 
are  based  on  program  code  generated  by  IBM  FORTRAN-Gl  and  FORTRAN-H 
compilers,  with  demonstration  runs  made  on  an  IBM  370/168  computer. 

Programs  SROW  and  DROV7  are  supplied  ready  to  execute  on  systems 
with  IBM-standard  FORTRAN  unit  members  for  the  card  reader  (unit  5) 
and  line  printer  (unit  6).  The  programs  may  be  modified  for  other 
systems  simply  by  changing  the  statements: 

KR  = 5 
KW  = 6 


to  the  proper  unit  numbers.  These  statements  appear  near  the 
beginning  of  the  SROW  and  DROIV  MAIN  programs. 

2 . 4 Model  Generation  and  Program  Flow 

The  panel  finite-element  models  are  subdivided  into  two  major 
components:  the  immediate  neighborhood  of  the  damaged  fastener 

hole  and  the  remainder  of  the  structure.  The  region  near  the 
damaged  hole  is  generated  by  procedure  RING  as  an  annulus  substruc- 
ture with  boundary  nodes  on  its  inner  and  outer  circumferences. 

The  remainder  of  the  structure  is  generated  by  procedure  SBL  or 
procedure  DBL  and  statically  condensed  for  coupling  to  the  outer 
circumference  of  the  ring. 
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2.4.1  Procedure  SBL 


Fig.  5 illustrates  the  hierarchy  and  numbering  convention  for 
the  substructure  generated  by  procedure  SBL.  Procedure  HOLEL  is 
used  to  generate  the  region  around  each  fastener  hole.  Two  calls 
to  HOLEL  are  required:  one  for  the  undamaged  holes  for  which  the 

inner  boundary  of  HOLEL  is  set  at  Rq  = RI , and  one  for  the  damaged 
hole  for  which  Rq  = 2.52  (RI)  to  accommodate  the  ring.  The  undamaged 
HOLEL  is  assembled  repeatedly.  Procedure  LUG  is  used  to  generate 
the  lower  portion  of  the  panel. 

The  following  node  numbering  scheme  has  been  adopted  to  facili- 
tate automatic  generation  of  the  substructure  interconnections  while 
satisfying  the  requirement  that  the  boundary  nodes  be  numbered  last 
to  permit  static  condensation: 

1.  The  HOLEL  inner  boundaries  are  numbered  first,  in 
the  lef t-to-right  sequence  corresponding  to  the 
hole  position.  Thus,  nodes  1 to  24  are  placed  in 
the  first  HOLEL,  nodes  25  to  48  in  the  second,  etc. 

2.  The  mid-edge  nodes  on  the  fastener  row  centerline, 
the  fastener  row  upper  edge  nodes,  and  the  fastener 
row  lower  edge*  nodes  are  numbered  in  order,  from 
right  to  left.  Fictitious  nodes  are  included  at 
the  hole  centers  to  facilitate  the  numbering 
algorithm.  These  are  later  restrained  without 
having  received  any  stiffnesses  during  the  assembly. 

3.  Finally,  24  additional  node  numbers  are  assigned 
to  the  inner  boundary  of  the  damaged  HOLEL,  and 
the  original  node  numbers  are  reassigned  to  the 
lower  edge  of  the  LUG  substructure.  Some  of  the 
latter  will  be  left  over,  and  are,  therefore, 
restrained  without  receiving  stiffnesses. 

The  node  numbers  shown  in  Fig.  5 are  the  final  assignments  for  a 

case  in  which  the  third  in  a ro w of  three  fastener  holes  is  damaged. 


The  lower  edge  node  numbers  are  also  assigned  to  the  upper  edge 
of  the  LUG  substructure . 
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The  various  substructures  are  now  assembled  and  the  fictitious 


nodes  are  restrained.  The  bottom  edge  of  the  panel  is  also  restrained 
and  nodal  forces  representing  the  bearing  loads  are  applied  at  the 
undamaged  holes,  as  indicated  in  Fig.  5.  Finally,  the  assembled 
substructure  is  statically  condensed  to  a "Cheshire  Cat"  (nothing 
left  but  the  smile) , consisting  of  the  24  nodes  on  the  damaged  HOLEL 
inner  boundary,  ready  for  final  assembly. 

2.4.2  Procedure  DBL 

Fig.  6 illustrates  the  hierarchy  and  numbering  convention  for 
the  substructure  generated  by  procedure  DBL.  The  model  creation 
is  generally  similar  to  what  has  just  been  described,  the  only 
significant  differences  being  some  details  of  the  node  numbering 
pattern  and  the  inclusion  of  QUAD4  elements  in  the  assembly  to 
account  for  the  stagger  between  the  two  fastener  rows.  The  end 
product  of  procedure  DBL  is  another  Cheshire  Cat,  ready  for  final 
assembly . 

2.4.3  Executive  Flow 

Programs  SROW  and  DR0V7  have  the  same  executive  flow,  summarized 
in  Fig.  7.  The  input  data  are  read  and  printed  for  reference,  and 
procedure  SBL  or  procedure  DBL  is  called  to  create  the  appropriate 
Cheshire  Cat.  The  Cheshire  Cat  is  designated  as  substructure  1 
for  the  final  assembly.  At  this  point,  the  process  of  setting  up 
the  final  substructure  is  started.  A re-entry  point  appears  in 
the  midst  of  the  interconnection  generation  algorithm  (statement 
205  in  either  program)  to  permit  the  band-margin  computation  sub- 
routine (ORK)  to  erase  the  K and  q data  from  the  previous  analysis. 
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The  Cheshire  Cat  is  now  assembled,  and  procedure  RING  is  used 
to  generate  the  cracked  ring  (substructure  2) , which  is  also  assembled. 
Nodal  forces  for  the  bearing  load  on  the  damaged  hole  are  now  applied, 
and  a global  displacement  solution  is  obtained  for  the  structure. 
Finally,  a back-substitution  solution  is  executed  for  the  cracked 
ring,  displacements  are  extracted  for  one  or  two  PCRK59  elements, 
and  NAS A/ASTM- standard  stress  intensity  factors  and  are  com- 

puted. Logical  IF  statements  at  the  end  of  the  loop  determine 
whether  the  angular  crack  position  is  to  be  incremented  for  another 
analysis.  If  another  case  remains  to  be  analyzed,  the  crack  angle 
is  incremented  and  the  program  branches  back  to  the  re-entry  point 
at  statement  205. 

2 . 5 Output  Conventions  and  Error  Messages 

Fig.  8 illustrates  a sample  output  from  program  DROW  with  all 
optional  printout  deleted.  The  heading  identifies  the  program 
which  was  executed  and  repeats  the  user's  input  data  for  checking. 

Below  the  heading  appears  a table  giving  the  angular  position, 
and  for  each  crack.  If  the  input  data  have  been  given  in  units 

of  lbs.  for  bearing  force,  psi  for  Young's  modulus  and  inches  for 
dimensions,  and  are  in  units  of  psi  /in.  The  crack  angles 

are  given  in  degrees.  The  output  from  program  SROW  is  similar. 

Error  messages  printed  by  procedure  RING  and  by  the  FEARL-2 
software  have  been  documented  elsewhere  [2,3].  In  addition,  pro- 
cedures SBL  and  DBL  print  error  messages  and  terminate  execution 
if  the  restriction  given  by  Eq.  5 is  violated.  Finally,  programs 
SROW  and  DROW  print  warning  messages,  change  the  panel  length,  and 
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continue  execution  if  the  user  has  input  a length  too  short  com- 
pared to  the  center-to-center  spacing  of  the  fasteners.  The  pro- 
grams are  always  executed  with: 


LENGTH  >_  2 (CLSPC)  for  program  SROW 
LENGTH  > 3 (CLSPC)  for  program  DROW 


(6) 


2 . 6 Status  of  Programs 

Programs  SROW  and  DROVJ  are  maintained  as  029-punched  FORTRAN-IV 
source  decks,  sequenced  by  individual  subroutine.  All  input  options 
described  in  subsections  2.1  and  2.2  have  been  exercised  successfully 
for  panels  containing  up  to  three  fastener  holes  per  row.  Source 
listings  appear  in  Appendix  A for  program  SROW  and  procedure  SBL, 
and  in  Appendix  B for  program  DROW  and  procedure  DBL . Listings 
for  the  other  procedures  appear  in  a previous  report  in  this 
series  [ 2 ] . 
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Section  3 


DEMONSTRATION  EXAMPLES 

Figure  9 illustrates  an  example  of  a single  row  of  three 
fastener  holes,  each  loaded  by  a 1,000-lb.  bearing  force  in  a 
1-inch  thick  panel.  Sample  analyses  of  this  structure  were  run 
with  the  center  and  right  holes  damaged  by  a single  crack.  Butter- 
fly plots  for  and  are  presented  in  Fig.  10  (center  hole 

damaged)  and  Fig.  11  (right  hole  damaged) . The  observed  behavior 
is  similar  to  solutions  for  cosine  bearing  loading  applied  to  an 
attachment  lug  detail  t 1 ] - In  particular,  maxima  for  are  seen 
to  occur  for  cracks  at  + 105  degrees  from  the  line  of  action  of  the 
applied  load  (9  = 195  and  345  degrees) , in  the  symmetric  case  (Fig.  10) , 
and  significant  values  occur  for  cracks  on  the  load  line  (9  = DO 
degrees)  due  to  the  local  hoop  tension  effect  of  the  bearing  pressure. 
Asymmetric  behavior  is  evident  for  both  and  K.  in  the  case  of 
the  right  hole  damaged  (Fig.  11) . A curious  result  for  this  case 
is  that  the  maximum  K ^ occurs  when  the  crack  is  positioned  away  from 
the  neighboring  fastener  hole  (9  = 345  degrees) . 

Figure  12  illustrates  an  example  problem,  for  a similar  struc- 
ture with  two  rows  of  fasteners.  The  results  of  two  example  analyses 
of  this  structure  are  plotted  in  Figs.  13  and  14.  In  Fig.  13,  the 
center  hole  in  the  upper  row  is  damaged.  These  results  are  comparable 
to  the  behavior  of  the  single-row  cases,  except  that  asymmetry  is 
now  introduced  by  the  presence  of  the  second  row  of  fasteners.  Also, 
the  maxima  of  now  shift  downward  toward  the  second  fastener  row. 

In  Fig.  14,  the  center  hole  in  the  lower  row  is  damaged.  Comparison 
with  Fig.  13  shows  that  and  have  increased  roughly  by  a factor 
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of  2,  an  effect  which  may  be  attributed  to  the  fact  that  the  damaged 
hole  is  located  deeper  in  the  load-transfer  zone  of  the  structure. 
Solution  asymmetry  is  very  slight  for  this  case. 
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Section  4 


DISCUSSION  AND  CONCLUSIONS 


Unfortunately,  there  exist  no  independent  solutions  for 
cracks  at  fastener  holes  loaded  in  bearing.  However,  some  approxi- 
mate correlations  can  be  made  by  comparison  with  solutions  available 
for  similar  problems,  if  the  bearing  forces  are  translated  into 
average  tension  for  the  nominal  cross  section  of  the  panel.  The 
conversions  are  given  by: 

o,  * 3 x 1,000  lb./30  in2  = 100  psi 

l (7 ) 

c?2  - 6 x 1,000  lb./35  in2  = 171.5  psi 
for  the  single-  and  double-row  cases,  respectively.  Thus,  one 
might  logically  expect  to  find,  e.g.,: 

KI  (Fig.  14 ) /KI  (Fig.  10)  = 1.715  (8) 

for  a given  crack  angle,  together  with  a similar  result  for  Kjj* 

The 'following  summary  of  calculations  tends  to  confirm  the  expected 
behavior : 


Item 

From  Fig.  10 

From  Fig.  14 

Ratio 

% Err 

KI 

at 

6 = 

0° 

250  psi  /In 

431  psi  /in 

1.72 

0.3 

KI 

at 

6 = 

180° 

250 

436 

1.74 

1.5 

II 

at 

e = 

210° 

181 

324 

1.79 

4.4 

II 

at 

e = 

330° 

181 

320 

1.77 

3.2 

An  approximate  comparison  may  also  be  made  with  published  solutions 
for  cracks  in  an  infinite  plate  under  uniform  tension  [4] : 


Kj  (6=0°) 


a/4b  sin  (t^) 


(9) 
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for  a row  of  repeated  cracks,  or: 

KI(0=O°)  = a/ Tra  F (a/Rj)  (10) 

for  a single  crack  emanating  from  a single  hole.  The  quantities 
a,b,c,RT  are  defined  in  Fig.  15.  For  the  present  case,  the  results 

as  follows: 

Computed  (Eq.  9)  (Eq.  10) 

250  246  271 

431-436  421  465 

where  the  stresses  a^,  given  by  Eqs.  7 have  been  used  in  Eq . 9 
and  Eq.  10.  The  fact  that  the  computed  values  are  below  the 
classical  solution  for  a crack  at  a hole  may  be  attributed  to 
incomplete  development  of  the  nominal  tension  in  the  immediate 
vicinity  of  the  fastener  bearing  load.  Also,  it  is  comforting  to 
note  that  the  f inite-element  results  are  bounded  by  the  two  classical 
solutions,  within  a range  of  about  10  percent. 

The  results  presented  in  this  report  demonstrate  that  the 
assumed-stress  hybrid  finite  element  method  is  capable  of  providing 
stress  intensity  factor  solutions  for  a wide  variety  of  structural 
details.  The  specific  program  codes  documented  in  this  report  treat 
panels  loaded  by  bearing  through  a single  row  or  a double  staggered 
row  of  fastener  holes  near  one  edge  of  the  panel.  The  programs  are 
presently  restricted  to  fairly  widely  spaced  fasteners,  in  order  to 
remain  in  the  parameter  ranges  for  which  the  most  accurate  performance 
of  the  individual  finite  elements  was  demonstrated  in  earlier  work. 


of  the  comparison  are 

Program 
SROW  (Fig.  10) 
DROW  (Fig.  14) 
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This  restriction  and  an  approximate  comparison  of  the  computed 
results  with  classical  solutions  indicate  that  the  finite-element 
computations  for  and  are  within  5 to  10  percent  of  the 

exact  answers.  A better  error  assessment  does  not  seem  likely  in 
the  near  future  because  of  the  difficulty  in  obtaining  classical 
solutions  for  the  very  complicated  geometry  and  loading  associated 
with  these  structural  details.  Only  a few  demonstration  examples 
have  been  presented,  but  the  possibilities  for  these  programs  are 
legion.  Situations  in  which  the  bearing  load  varies  from  one 
fastener  to  another  may  be  of  interest  to  the  detail  designer.  One 
extremely  interesting  example  is  the  question  of  stress  intensity 
at  a hole  which  has  cracked  sufficiently  to  unload  its  fastener. 
Will  this  result  in  crack  arrest,  or  will  the  bearing  loads  from 
nearby  fasteners  still  impose  a K great  enough  to  continue  crack 
growth  or  to  cause  fracture?  This  and  other  example  solutions  will 
be  presented  in  the  final  report  in  this  series. 
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FIG.  5 SUBSTRUCTURE  HIERARCHY  (PROCEDURE  SBL) 
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?IG.  6 SUBSTRUCTURE  HIERARCHY  (PROCEDURE  DBL) 
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SAMPLE  OUTPUT  (PROGRAM  DROW) 


— -A Kit  I 1 1 1 

0 .1  .2  .3 

KSl/ltf 


FIG.  10  BUTTERFLY  PLOT  (CENTER  HOLE  DAMAGED) 
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FIG.  12  EXAMPLE  OF  PANEL  WITH  TWO  ROWS  OF  THREE 
FASTENER  HOLES 
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FIG.  15  DEFINITION  OF  CRACK  PARAMETERS  FOR  CLASSICAL  SOLUTIONS 
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PRINT  OUT  PROGRAM  STATISTICS  TO  DATE:  SROW0109 

IF (KT1  .EQ.  KW  ) GO  TO  65  SROWOllO 

WRITE  (KW»6007)  WIDTH*  SM(.l*ETA  SROWOlll 
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